A simple model for dynamical melting of moving vortex 
lattices interacting with periodic pinning 
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Melting of moving vortex lattices in clean superconducting films with periodic pinning is studied 
by a dynamical "cage" model, based on a mean-field treatment of Langevin's equations for the whole 
vortex lattice, assuming elastic fiow. In the frame moving with the velocity of the vortex lattice 
center of mass, the model describes forced vibrations of a single vortex tied to a spring. The vortex 
displacements due to thermal fiuctuations and to the periodic pinning force, and the relationship be- 
tween the driving force and the vortex velocity (V-I curves) are obtained by a perturbation method, 
valid for high velocities and for both weak and strong periodic pinning. The dynamical melting 
temperature is calculated as a function of the vortex velocity using Lindemann's criterion. Appli- 
cation to a square defect lattice gives, for strong pinning, dynamical melting lines and anisotropic 
V-I curves that agree qualitatively with numerical results. 
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I. INTRODUCTION 



The study of dynamical phases of lattice structures 
moving in a periodic medium has received a great deal of 
attention lately due to its relevance to physical systems 
such as adsorbed atomic layers in boundary JubricationEJ, 
vortices in Josephson Junction Arrays ( JJA)0 and in films 
with artificial defect latticed^. 

Numerical studies of driven vortices in clean supercon- 
ducting films intcmcting with a.neriodic lattice of colum- 
nar defects (CD)cnJ and in JJAtlQ find dynamical melting 
of a moving vortex-lattice (VL) into a moving vortex liq- 
uid. In these systems the moving vortices order in a VL 
at low temperatures when the driving force magnitude 
is sufficiently large. Melting of the VL is found to take 
place at a temperature that, for a given vortex density, 
decreases with the driving force magnitude, approaching 
the equilibrium VL melting temperature when this mag- 
nitude approaches infinity. These studies also show that 
the vortices center of mass (CM) velocity is in general not 
parallel to the driving force. This leads to anisotropic V- 
I curves, with voltages both longitudinal and transverse 
to the direction of the applied current. 

In this paper we show that dynamical melting lines 
and anisotropic V-I curves similar to those found in the 
above mentioned numerical studies can be obtained using 
a simple model for the vortex dynamics. 

In driven vortex systems interacting with random or 
periodic pinning, moving vortex lattices are expected to 
occur at low temperatures and in the limit of very large 
CM velocities, because the pinning potential is averaged 
in the direction of motion, as first pointed out by Schmid 
and Haugem. For periodic pinning, Ref. shows that in 



this limit the moving vortex lattice can be commensurate 
or incommensurate with the periodic pinning potential, 
depending on the direction of motion and on the vortex 
density. Commensurate lattices occur only for motion 
along directions of high symmetry of the pinning poten- 
tial, such as [1,0], [0,1], [1,1] and [-1,1] for scuiap: sym- 
metry, if the vortex density is not too highoiJ'E3. For 
motion along other directions, the moving vortex lattices 
are incommensurate. Recent numerical simulations find 
that both moving commensurate and incommensurate 
lattices .jp^y undergo dynamical melting into a moving 
liquidLTLrEj. However, dynamical melting of commensu- 
rate lattices is much more complicated than that of the 
incommensurate ones. The reason is that the moving 
commensurate lattices are pinned in channels running 
along the direction of motion (transverse pinning), and 
melting into a liquid ocf|U|Ca simultaneously with depin- 
ning from the channelsQ&tj. In this paper we consider 
only dynamical melting of incommensurate lattices. 

Dynamical melting of vortex lattices is still poorly un- 
derstood. It is unclear from the simulation results re- 
ported so far if the change from lattice to liquid spa- 
tial symmetries is a true phase transition, with accompa- 
nying singularities in correlation functions, or merely a 
crossover. In this paper we do not offer any insight into 
this fundamental question. Our aim is to estimate the dy- 
namic melting temperature based on a well know method 
for the equilibrium melting temperature, namely Linde- 
mann's criterion. This method was proposedin a classi- 
cal paper by Houghton, Pelcovitz and Sudb0Eil. They es- 
timate the equilibrium VL melting temperature for clean 
bulk superconductors by applying Lindemann's criterion 
to the mean-square vortex-line displacement caused by 
thermal fiuctuations. A simple model to estimate the 



1 



same temnerature was proposed later by Frey, Nelson 
and FishertJ. It considers a single vortex line trapped 
in a harmonic potential representing its interaction with 
the other vortex lines in the lattice, assumed straight 
and fixed at their equilibrium lattice positions. We refer 
to this model here as the equilibrium cage model. The 
equilibrium melting temperature is estimated by apply- 
ing Lindemann's criterion to the mean-square thermal 
fluctuations of the caged vortex line. The result agrees 
with that of Ref. |ll|, if the harmonic potential is properly 
choscnO. This simple idea can also be applied to two- 
dimensional VL melting in thin films. In this case the 
equilibrium cage model consists of a point vortex tied 
to a spring, representing its interaction with the other 
point vortices, assumed fixed at their lattice positions. 
The melting temperature is obtained by applying Linde- 
mann's criterion to the caged vortex mean square dis- 
placement due to thermal fluctuationsEJ. The resulting 
melting temperature agrees isath the one obtained from 
dislocation unbinding theoryc^, if the spring constant is 
properly chosen. Our model, which we call the dynamical 
cage model, generalizes the two-dimensional equilibrium 
one for a moving vortex lattice. The dynamical melting 
temperature is estimated by applying Lindemann's crite- 
rion to the mean-square vortex displacements caused by 
thermal fluctuations and by motion in the periodic pin- 
ning potential. Application of Lindemann's criterion to 
estimate the dynamical melting line of vortices in bulk su- 
perconductors interacting with randmn pinning was first 
proposed by Koshelev and Vinokuro 

We start from Langevin's equations describing two- 
dimensional vortices at temperature T, interacting be- 
tween themselves and with a periodic pinning potential, 
and driven by an external force. We consider vortex mo- 
tion at large CM velocities and make the following as- 
sumptions: i) The moving vortices order in a lattice, in- 
commensurate with the periodic pinning potential, ii) 
The vortex lattice flows elastically. This means that 
there are no dislocations and that vortex-vortex inter- 
actions can be approximated by harmonic forces. We 
then propose a mean-fleld like treatment of the resulting 
equations of motion. This approximation reduces the 
problem to that of a single vortex tied to a spring, both 
moving with the CM velocity, and interacting with the 
periodic pinning potential. We find that the dynamical 
cage model is governed by two equations of motion. One 
for the vortex displacement in the frame moving with the 
vortex lattice CM velocity (CM frame), and another re- 
lating the driving force and the CM velocity. The vortex 
displacement equation is found to be nonlinear, unlike 
that for the equilibrium cage model, due to the vortex- 
pinning interaction. These equations are solved by a per- 
turbation nifithod, similar to that introduced by Schmid 
and HaugerB for weak random pinning forces. It consists 
in expanding the time-dependent periodic pinning force 
(in the CM frame) in powers of the vortex displacement. 
We find that the perturbation expansion is valid at large 
CM velocities for both weak and strong pinning. From 



these solutions we obtain expressions for the vortex mean 
square displacement and for the relationship between the 
CM velocity and the driving force. The dynamical melt- 
ing temperature is obtained by applying Lindemann's 
criterion to the caged vortex mean-square displacement. 
These results are applied to a square pinning array, and 
compared to numerical simulation ones. 

This paper is organized as follows. In Sec. |l| we derive 
the dynamical cage model equations of motion and their 
perturbation theory solution. In Sec. III. we apply these 
results to strong pinning by a square array of pinning 
centers. The conclusions of the paper are presented in 



Sec. IV. Details of the perturbation calculation are given 



in the Appendix. 



II. DYNAMICAL CAGE MODEL 

We consider Ny two-dimensional vortices at tempera- 
ture T interacting between themselves and with a peri- 
odic potential produced by an array of pinning centers, 
and driven by a force id- The equations of motion arelij 



'7^ =fd + FJ-^ + F-f(r,) + r, , 
where j — 1, Ny, rj is the friction coefficient. 



(1) 



(2) 



is the force of interaction with other vortices, U^~^{r) 
being the vortex-vortex interaction potential in two di- 
mensions, 



F''-^(r,) = EHQ)C/QVQ- 
Q 



(3) 



is the force of interaction with the periodic pinning array, 
Q denotes the pinning array reciprocal lattice vectors, 
and Uq^ is the Fourier transform of the vortex-single 
pin interaction potential, and Tj is the random force, 
satisfying 



(r,a(t)r,;3(t')) = 2kBTi^5,^i5^,p5{t - t') 



(4) 



where (• • •) denotes average over the random force distri- 
bution. 

We assume that the vortices center of mass (CM) 
moves with constant velocity v, defined, as usual, by 



Ny ^ dt 



(5) 



Using Eq. (m) , we find that 
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(6) 
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To derive Eq. (g) we use the fact that for large N^, the 
random force term is neghgible, since ^ 



According to Eq. (|6|) , in order that v is constant must 
depend on time as weU as on the random force. To ob- 
tain physical results we average Eq. (|^) over time and 
over the random force distribution, that is 



1 1 r 



(7) 



where F^, a constant vector, denotes the average of 
over time and over the random force distribution and r 
is a time large compared with the characteristic times in 
Eq. (|l|). We interpret as the force due to the appHed 
current, so that Eqs. (j^) give the V-I curves. 

Now we consider a vortex-lattice flowing elastically and 
transform to the CM frame. We write the vortex posi- 
tions as 



(8) 



where Rj and Uj{t) {j — l,...,Ny) are, respectively, 
the vortex equilibrium positions and displacements from 
equilibrium in the CM frame. In elastic flow the displace- 
ments Uj (t) are small and the vortex- vortex interactions 
can be approximated by harmonic forces. Our dynamical 
cage model considers a particular vortex, and substitute 
its elastic interaction with the other vortices in the lat- 
tice by a spring of constant k. The equation of motion 
for this vortex in the CM frame is found to be 

dui{t) 

v—jT- = -^Mt) 



dt 



fd - yyv + r, . (9) 



The relationship between the spring constant and the 
vortex-vortex interaction potential will be presented 
shortly. Our dynamical cage model approachis similar to 
the mean-field theory introduced by Fishcrll3 for sliding 
charge-density waves interacting with random pinning. 

Our objective here is to use Eq. (^, together with Eq. 

, to calculate the vortex mean-square displacement re- 
quired for Lindemann's criterion, and to obtain the v vs. 
Fd relationship (V-I curves). First we note that Eq. (|^) 
depends on the particular vortex I, through the F^"^- 
term. Since F^'^ is periodic in the pinning-array lattice, 
its argument in Eq. (^) can be reduced, at any instant 
t, to a position within the pinning-array primitive unit 
cell. As mentioned in Sec. ^, we consider only moving 
incommensurate vortex lattices. In this case the reduced 
vortex positions are uniformly distributed over the unit 
cell area. Accordingly, we define the vortex mean square 
displacement as 



1 1 r 



(10) 



According to the above considerations, our dynamical 
cage model is described by the set of equations 



dii^ = -«uKi) + F-f (R/ + ui{t) + vi) 
dt 



^5^F''-''(R,+u,(t)+vt)+r,, (11) 



and 



''^ = + 1^ E 7 (F''-^(R, + u^t)) . (12) 



To obtain ui{t) we solve Eq. ( |ll| ) by the perturbation 
theory method introduced in Ref. |^, which consists in 
expanding F""'' in powers of Ui{t) and solving the result- 
ing equations by iteration. We keep only terms to first 
order. The relationship between v and F^ are obtained 
from Eq. (|l2|) by a similar expansion. The details of these 
calculation are given in the Appendix. 

For the mean square displacement we obtain to leading 
order 



2kBT 



2kBT 



E 

Q 

E 

Q 



Q'\Uc 



77-^(Q • v)-^ -|- 



Q^\Uq\ 



vHQ- 



4^2 



(13) 



To the same order of approximation the v vs. F^ re- 
lationship is 



rj ^ 772(Q • v)2 + k2 
Q 



Q V 



K ^ 772(Q. v)2 -I- 4^2^ 

Q 



(14) 



According to Lindemann's criterion the VL melts when 
= cj^a^, where a„ is the VL lattice parameter. Thus, 
our dynamical cage model only makes sense as long as 
< c|a2. This requires that in Eq. ( p^ ) both /cbT/k, 
and the velocity dependent terms are sufficiently small 
compared to a^. The latter are so for large enough v. 
To state more precisely the conditions under which the 
perturbation expansion in Eq. ( p^ is valid, we restrict 
our considerations to periodic pinning forces that can 
be represented by Eq. (|^) with only a few non-zero Uq. 
Typically those for Q corresponding to nearest and next- 
nearest neighbor reciprocal lattice points. In these cases 
Q ~ 27r/ap, ap being the pinning array lattice parameter. 
We also assume that the number of vortices is comparable 
to the number of pins, so that ^ Up. Under these as- 
sumptions, the condition for the smallness of both veloc- 



ity dependent terms in Eq. (13) depends on the relative 
strengths of the pinning and elastic forces. For weak pin- 
ning, defined as \Uq\ ^ Kap/2TT, these terms are always 
small, independent of v, so that the effects of vortex 
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FIG. 1. (a) Pinning lattice, coordinate system and angles 
defining the directions of the center of mass velocity, v, and 
the driving force, Fd- (b) Square pinning potential (Eq. (p^) 
with U2 = Ui/2. 

motion on are negligible. For strong pinning, that 
is \Uo\ ~ Kap/27r, the velocity dependent terms in Eq. 



( |13| ) are small for all Q, if t^wcosiJq ^ nap /2tt , where 
i9q is the angle between the particular Q and v. Sim- 
ilar considerations apply to the perturbation expansion 
for V, Eq. (jlj). For strong pinning, the perturbation ex- 
pansion breaks down if, for some Q entering the sums 
in Eq. (|l3|), Q • v = (cosi9q =0). In this case the 
velocity dependent terms in Eq. ( ]l3| ) are not small. This 
breakdown can also be seen from Eq. (^). If Q • v = the 
pining force F''"P(R/ -|-u/(t) -I- vt) has a static component, 
periodic in the direction perpendicular to v, which gives 
rise to a static vortex displacement comparable to for 
strong pinning. Such a displacement indicates that our 
assumption that the moving VL is incommensurate with 
the pinning array lattice breaks down for motion along 
the directions for which Q • v = 0. There are only a few 
directions of motion for which Q • v = since, as men- 
tioned above, only a few Q enter the sums in Eq. (13). 
We note that the terms in sum in Eq. ( p^ with Q • v 
do not contribute to v, so that there is no breakdown in 
the perturbation expansion for this quantity. 

In the limit of very large v {v 00) our dynamical cage 
model reduces to the static one, since vortex motion av- 
erages out the pinning potential (we assume Q • v 7^ 0) . 
In this limit it^ reduces to the equilibrium mean square 



displacement of a vortex tied to a spring. 



,2 _ 2kBT 



Ap- 



plication of Lindemann's criterion gives for the melting 
temperature ksTm = na^c^/i. To compare this esti- 
mate of Tm with the dislocation unbinding theory one, 
n is identified with k = (9^t/^^"(r)/9r^)r-=a„, where 
U^-"{r) = (0g/87rAIln(r/C) is the 2D vortex- vortex 
interaction potentiaO (00 is the flux quantum and A 
the film effective penetration depth ). The result is 



ory onelii if ~ 0.12. From here on we adopt for k and 
cl the above mentioned values. 

Application of Lindemann's criterion to Eq. ( [l3| ) gives 
the velocity dependent melting temperature Tdm{"^) as 



1 + Q^IUqI^ 



where 



kBTppiv)^^J2^ 



(15) 



(16) 



is identified with an effective temperature, resulting from 
vortex vibrations due to the periodic pinning force, that 
adds to the thermodynamic one in the expression for u^. 
The temperature Tpp is similar to the 'shaking tempera- 
ture' introduced by Koshelev and Vinokurfij for random 
pinning. It arises from a T-independent contribution to 
u^, the second term in the right hand side of Eq. (p^). 
There is also a T-dependent contribution to u'^ in Eq. 
( [Tsl ) that cannot be identified with an effective temper- 
ature. However, this term leads to a small correction 
in Tdmi'y), since the denominator in Eq. (|l^) is close to 
one. Thus, the dynamical melting condition is essentially 
that kB{Tdm{v) + Tpp{v)) = Ka^c|/2. The v dependence 
of rdm(v), Eq. (p^, is such that Tdm {"v) decreases with 
decreasing v. 

The V-I curves follow from the relationship between v 
and Fd, Eq. (p^. In the v ^ 00 Hmit v = Fd/rj. For 
finite v, vortex motion is no longer along the drive di- 
rection, leading to anisotropic V-I curves. The v vs. Fd 
relationship also allows us to obtain the melting temper- 
ature as a function of the driving force, which is more 
usual. 

In the next section we study in detail the dynamical 
cage model predictions by applying them to a typical pe- 
riodic pinning potential and carrying out numerical cal- 
culations. 



III. SQUARE PINNING-LATTICE 

We consider a square pinning array with lattice param- 
eter Up, and assume that the vortex-pinning interaction 
force is given by Eq. (^ with the Uq^ chosen as 



Ui for Q = ±^x, ±^y; 
U2 ior Q = ±^{x±y); 
otherwise. 



(17) 



The X and y axis are along the pinning lattice [1, 0], [0, 1] 
directions, respectively, as shown in Fig. |l|(a). We choose 
U2 = Ui/2. This gives the pinning potential shown in 
Fig. |l](b). We set Ui/kuI = 11211 = 0.16. As discussed 
temperature agrees with the dislocation unbinding the- in Sec. |l[ this value corresponds to strong pinning. 



K = 0Q/87rAa„ and ksTm = CL4'tl'^^''^J^- This melting 
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FIG. 2. Center of mass velocity components, Vx and Vy, 
and root-mean-square caged vortex fluctuations, u, vs. driv- 
ing force magnitude, Fd, for several driving force directions a 
at T = 0. (a) 32.5° < a < 40°: vortex motion is trapped in 
the [1,1] direction {vx = Vy) at the Fd values indicated by the 
vertical dot-dashed lines in the top panel, (b) < a < 32.5: 
vortex motion is trapped in the [1,0] direction {vy — 0). 
The horizontal lines in the bottom panels of (a) and (b) are 
u = clOi, for B = 2_B^ and cl = 0.12 (see text), where the 
VL melts. 



We choose to present the numerical results for the ex- 
pressions derived in Sec. || as functions of the driving 
force Fd: instead of v. To do so we invert Eq. ( p^ by 
an iterative method. We start at a high value of the 
driving force magnitude, Fd = bOnap, and a given ori- 
entation of Fd vifith respect to the x-axis, a (Fig. 0(a)), 
and approximate the solution by the result for infinite 
drive v — Fd/ij. Then Fd is decreased by small steps, 
keeping a fixed. At each step, v in the right hand side 
of Eq. ( p^ ) is replaced by the v obtained in the previ- 
ous step. From the v vs. relationship thus obtained, 
the u and Tdm vs. Fd curves follow using Eqs. (|l3|) and 
(p5|). The results are shown in Figs. In the fig- 

ures shown here the following units are used, i) Force: 
Fmax — Snlli/ap — magnitude of the maximum force 



exerted by the square pinning potential defined in Eqs. 
(|l7|). ii) Velocities: Fmax/v- iii) Temperatures: = 
VL equilibrium melting temperature (see Sec. 

In Fig. 1^ we show v and u as functions of Fd for sev- 
eral a at T = 0. The horizontal lines in the bottom 
panel of Figs. |^(a) and |^(b) indicate the u — CLiy lines 
for a vortex density corresponding to two vortices per pin 
{B — 25^), with Cl — 0.12, at which the VL melts, ac- 
cording to Lindemann's criterion. The region of validity 
of our dynamical cage model is below this line. The v 
vs. Fd curves in the top panels of Figs. |2|(a) and ^(b) 
show that the CM velocity is not parallel to the driv- 
ing force, which leads to anisotropic V-I curves. This 
anisotropy is more clearly seen by plotting the direction 
of motion, 9, as a function of the driving force magni- 
tude, as shown in Fig. |[ for fixed directions of drive, 
a. These plots show that, depending on a, the mov- 
ing VL is attracted towards the [1,0] or [1,1] directions. 
This behavior is similar to that seen in recent numerical 
simulationaJ. Above the cl = 0.12 lines in the bottom 
panel of Figs. |^(a) and^(b) we find a discontinuous jump 
in u, accompanied by the trapping of moving caged vor- 
tex along [1,0] for < a < 32.5° [Fig. |(b)], or along [1,1] 
for 32.5° < a < 45° [Fig. ||(a)]. This trapping is similar 
to transverse pinning observed in numerical simulations 
and JJA experiments.u Q This is interesting because it 
indicates that, even though our dynamical cage model 
is not strictly valid in the region where trapping occurs, 
it contains physical ingredients capable of describing the 
phenomenon. We also note that the jump in u occurs 
because for vortex motion along [1,0] and [1,1] there are 
terms in the sums in Eq. ([ij) for which Q • v = 0. Ac- 
cording to the discussion in Sec. |l| the perturbation ex- 
pansion for u breaks down and a large contribution to 
u results from the terms in the sum in Eq. ( p^ ) with 
Qv^O. 
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FIG. 3. Direction of motion, 9, as a function of the driving 
force magnitude for drive directions in the range 5° < a < 40° 
in regular increments of 5°. The arrow indicates increasing a. 
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FIG. 4. Temperature dependence of: (a) center of mass 
velocity y-component and (b) root-mean- square displacement 



curves, u vs. 



Fd. 



FIG. 5. Top: dynamical melting lines for a in the range 
10° < a < 40°. Bottom: dynamical phase diagram Fd vs. a 
for T = . 



The finite-T behavior of the results is shown in Fig. 
^. We find that the v vs. Fd curves change little with 
T as illustrated in Fig. ^(a). We also find that in 
the T-dependent term in Eq. ( p^ ) gives a small contri- 
bution. This is so because the T-dependence in u and 
V results from the third terms in the right hand side of 
Eqs. disl ) and (p^, which are, essentially, the product of 
two quantities which are small in the region of validity of 
our dynamical cage model. 

The dynamical melting lines Tdm vs. Fd are shown in 
the top panel of Fig. |^ for several a. In the bottom panel, 
the dynamical phase diagram Fd vs a for T = is shown. 
The Fd vs. a dynamical melting line in the bottom panel 
corresponds to the intersections of the Tdm vs. Fd curves 
shown in Fig. || top with the T — line. We find that for 
other temperatures T < Tm the Fd vs a phase diagrams 
are similar. 



IV. CONCLUSION 

In conclusion then, we introduce a simple dynamical 
model for vortex lattices flowing elastically in the pres- 
ence of periodic pinning. The model is solved by a per- 
turbation method that allows the calculation of the dy- 
namical melting curves, using Lindemann's criterion. 



and of the V-I curves. This approach is valid for peri- 
odic pinning potentials in general, and predicts that the 
moving vortex-lattice melting temperature varies signif- 
icantly with the CM velocity only for strong pinning. 
The perturbation solution is applied to a square pinning 
lattice. The results shovi. mialitative agreement with nu- 
merical simulation onescTD. Namely: i) the predicted 
dynamical melting temperature approaches the thermal 
equilibrium one in the limit of very large center of mass 
velocities, and decreases with the velocity. The model 
predicts that this decrease is, essentially, as ii) The 
center of mass velocity is not parallel to the driving force 
in general, leading to anisotropic V-I curves. The vortex 
lattice is attracted towards the [1,0] or [1,1] direction, 
depending on the direction of drive. Transverse pinning, 
with the vortex lattice motion restricted to the [1,0] or 
[1,1] directions, is found outside the region where the 
model is strictly valid. 

For other periodic pinning potentials, and for strong 
pinning, we expect similar dynamical melting lines and 
anisotropic V-I curves. 

It is significant that the model predicts anisotropic V- 
I cuH/es so closely resembling the numerical simulation 
onesd. These curves follow from the relationship between 
the CM velocity and the driving force, Eq. (1^). The 
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only approximation used in our solution of this equation 
is the mean-field treatment of Langevin's equations for 
the whole vortex lattice. Our results for the V-I curves 
show that this approximation is justifiable. 

The model can be also applied along the lines described 
here to study the effects of weak random potentials, ei- 
ther from material defects or from imperfections in the 
pinning array, on the melting lines and V-I curves. It 
can also be generalized to study transverse pinning more 
rigorously. Work along these lines is under way and will 
be reported elsewhere. 
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APPENDIX A: 

To obtain the perturbation expansion of Eq. ( |ll|) we 
write it, using Eq. (||), as 



duijt) 
' dt 



Q 



- Y,(NvSj,i - i)e*Q-[R.+".(*)l + Ti , (Al) 



where Fq'' = — iQC/q^ and luq = Q • v. Expanding 
giQ uj £j,g^ order, and carrying out the Pourier trans- 
formation in time, assuming periodic boundary condi- 
tions in the time interval (0, r), justifiable for large r 
(r ^ CLp/v), we find 



-irjLU + K _ 
W J — J- 



where u^(a;) is the thermal displacement 



(A2) 



(A3) 



P;(w) is the Pourier transform of the random force and 
6k denotes the Kronecker delta. The first order solution 
f or u ijuj) is obtained by neglecting the linear term in Eq. 
@. We find 



u«M=urH-H5: 



Q 



ir]iUQ + K 



(A4) 



To derive Eq. (|A4|) we use e'^'^^ = (5q,o, since the 



vortex and pinning lattices are incommensurate. 



The second order solution is obtained by substituting 



Uj in the linear term in the left hand side of Eq. ( A2) by 
the first order solution. The result is 



Q 



— J2{NJj,i - l)e''^-^HQ ■ uJ{lj + ljq) 



E 

Q,Q' 



-p^-P^rt ■i7if-Pr„iK-R, r 

{—irjULi + K){i?]ujQf + k) 



(A5) 



where ujq> = Q' • v, ljk = K • v, and K = Q + Q'. With 
this solution, one can calculate the mean square displace- 
ment for a particular vortex I, uf = ^ /J"(|ui(t)p), which 
up to the second order in the pinning potential is ex- 
pressed by: 



2kBT 



E 

Q,Q' 



T^v-p -pu-p iK-R, 



2ksT 1 g%-^-F-f 

[(iV, -2)e'K-R. +jq _q,] . 



(A6) 



Taking the average over the Ny vortices and assuming 
iV„ > 1, Eq. ^ is derived. 

A similar procedure is used to calculate the driving 
force dependence of v. Expanding Eq. in the form 



£ E (Fq'e*'^-^^- e^-«* [l + Q • u, (t)] ) , (A7) 



and using the second order solution for ui{t) we obtain 
Eq. (0). 
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